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We present a summary of results of the joint CP-PACS and JLQCD project toward a 2+1 flavor full 
QCD simulation with the (a) -improved Wilson quark formalism and the Iwasaki gauge action. 
Configurations were generated during 2002-2005 at three lattice spacings, a ~ 0.076, 0.100 and 
0.122 fm, keeping the physical volume constant at (2.0 fm) 3 . Up and down quark masses are 
taken in the range mps/my ~ 0.6 — 0.78. We have completed the analysis for the light meson 
spectrum and quark masses in the continuum limit using the full configuration set. The predicted 
meson masses reproduce experimental values in the continuum limit at a 1 % level. The average 
up and down, and strange quark masses turn out to be m^ s (/i = 2 GeV) = 3 .50(14) (ijf) MeV 
and m^ s (ii = 2 GeV) = 9 1.8 (3. 9) (tfl) MeV. We discuss our future strategy toward definitive 
results on hadron spectroscopy with the Wilson-clover formalism. 
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1. Introduction 

The calculation of the light hadron spectrum and quark masses is a fundamental and necessary 
step for the entire area of lattice QCD simulations. Moving from a precision quenched calcula- 
tion [@ to two flavor (N f = 2) full QCD [§], we observed a significant dynamical up and down 
("light") quark effect, which removes most of the 0(10%) systematic deviation in the quenched 
QCD spectrum from experiment. It was also found that light quark masses are significantly re- 
duced in Nf = 2 full QCD. In order to remove the quenching error of the heavier strange quark, 
the CP-PACS and JLQCD collaborations have jointly pursued a 2+1 flavor (Nf = 2+ 1) full QCD 
simulation [[J since 2001. We employ the Wilson quark formalism as in our quenched and Nf = 2 
studies, preferring an unambiguous quark-flavor interpretation over the computational ease of the 
staggered formalism carried out by the MILC collaboration [Q]. 

The project explored the light quark mass range corresponding to m n jm p = 0.6 — 0.78 for 
which the configuration generation has been completed at three lattice spacings in the fall of 2005. 
In this article, we present a summary of the results for light meson masses, light quark masses, 
pseudoscalar (PS) decay constants and the Sommer scale, evaluated in the continuum limit using 
the full set of configurations. We also discuss systematic error from chiral extrapolations by com- 
paring fits with polynomial functions of quark masses and those based on chiral perturbation theory 
CfFT). 

2. Production of the gauge configuration 

For the lattice action, we employ the renormalization group (RG) improved Iwasaki gauge ac- 
tion and the clover quark action with the improvement coefficient c$w determined non-perturbatively 
for the RG action The choice of the gauge action is made to avoid a first-order phase transition 
(lattice artifact) observed for the plaquette gauge action [§J. 

Configurations are generated with the Polynomial Hybrid Monte Carlo (PHMC) algorithm. 
(See [Q] for our implementation.) The molecular dynamics time step St and the polynomial order 
Np iy are chosen such that the HMC and the global Metropolis acceptance rate achieves 85% and 
90%, respectively. 

Simulations are performed at three values of the coupling constant such that a 2 is placed at an 
even interval. The physical volume is fixed at (2.0 fm) 3 . Main simulation parameters are listed in 
Table [j]. At each coupling, we generate configurations for ten combinations of hopping parameters 
(k uc i,k s ), five for the ud quark mass taken in the range of the pseudoscalar (PS) to vector (V) 
meson mass ratio of nips/my ~ 0.6 — 0.78 and two for the strange quark mass chosen around 
nips/mv ~ 0.7. 





size 


a [fm] (i^-input) 


a [fm] (0 -input) 


trajectory 


1.83 


16 3 x 32 


0.1222(17) 


0.1233(20) 


7000-8600 


1.90 


20 3 x 40 


0.0994(19) 


0.0995(19) 


5000-9200 


2.05 


28 3 x 56 


0.0693(26) 


0.0695(26) 


6000-6500 



Table 1: Main simulation parameters. 
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3. Measurement and analysis 

Measurements are made at every 10 HMC trajectories. We use the combination of smeared 
source and point sink, with which the effective masses reach plateau earliest and the statistical 
errors are smallest. Meson masses and amplitude of the correlator are obtained from single mass 
X 2 fits to correlators (P(t)P(O)), (V(f)V(O)) and (A 4 (f )P(0)>, where P, V and denote the PS, the 
vector and the non-perturbatively 0(a)-improved [Js[] axial-vector current, respectively. We include 
correlations in time but ignore those among correlators, since our statistics are not sufficient for 
the latter. Errors are estimated by the binned jackknife method with the bin size of 100 HMC 
trajectories. 

For chiral extrapolations of meson masses, we use two definitions of the quark mass, the vector 
Ward identity (VWI) quark mass m q WI = (1/ff— l/fQ.)/2 (k^ is the critical hopping parameter 
where mps at k uc i = fc s = K va i = K c vanishes), and the axial-vector Ward identity (AWI) quark mass 

V^A^(x)/ (2P(x)). Quark masses are determined from chiral fits with m q , because m„ 
shows large scaling violation. 1 We use chiral fits with m q WI for other quantities. 

Chiral fits are made to light-light, light-strange and strange-strange meson masses simulta- 
neously ignoring their correlations, using a quadratic polynomial function of the sea and valence 
quark masses; 



where Ms = diag(m„ ( ;,m (1 j,m i ), My = diag(m m n,m va i2), and "tr" means trace of matrices. We set 
A = for fits of m 2 PS with m v q WI , while A=B S = C S \= C S2 = for those with mf 1 . 

The physical point is fixed for two cases. The "^T-input" takes the experimental values m % = 
0.1350 GeV, m p = 0.7684 GeV and m K = 0.4977 GeV as inputs. In the "0-input" case, m n , m p 
and m$ = 1.0194 GeV are taken as inputs. The lattice spacings determined from the K- and <j>- 
inputs are consistent with each other, as shown in Table fiL 

Static potential V (r) is determined from Wilson loops for smeared gauge links and is fitted to 
a form V(r) = c — a/r + or. The Sommer scale ro calculated from the fit is extrapolated to the 
physical point as 1/ro = A + BstrM s using ni^ WI . 

4. Physics results 

4.1 Light meson spectrum 

The meson masses are well fitted by linear functions in a 2 , as shown in Fig. [j]. We obtain in 
the continuum limit, 



The VWI quark mass for the ud quarks is negative at our simulation points. This originates from a lack of chiral 
symmetry of the Wilson quark action. This is another reason to prefer the AWI definition. 



/(M S ,M V ) = A + B s trM s + Bi/trM v + D sy trM s trM v 

+C S itrM2 +C 52 (trM s ) 2 + CyitrM2 + C^ 2 (trM v ) 2 



(3.1) 



m K * =0.8961(72) GeV, m = 1.023(14) GeV input), 
m K = 0.495(14) GeV, m K * = 0.8947(12) GeV (0 - input). 



(4.1) 
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Figure 1: Continuum extrapolation of meson masses, comparing with those for quenched and Nf = 2 
QCD t^]. Note that the quenched and Nf = 2 simulations are made with the one-loop perturbatively 0(a)- 
improved clover action. Thus extrapolations are made linearly in a. 



These values are consistent with experiment within the quoted statistical errors of about 1% level. 
The errors are larger than those achieved in the Nf = 2 case, and possible deviations from Nf = 2 
QCD are not detected. We observe that the scaling violation in Nf = 2 + 1 QCD seems to be larger 
than that in the quenched or Nf = 2 case [|2|] even though the actions used in these cases are less 
improved than the present Nf = 2 + 1 case. We may argue, however, that the magnitude of the 
scaling violation in Nf = 2 + 1 QCD is not particularly large. A quadratic fit of the K* meson 
mass with /f-input mg* = «?o (1 +c(A QCD -a) 2 ) yields m = 8960(73) MeV and c = -1.29(59) for 
Aqcd = 200 MeV. The 0(1) magnitude of the coefficient c is reasonable. 



4.2 Quark masses 

The physical quark mass is determined for the MS scheme at the scale \i = 2 GeV. Lattice 
results are translated to the MS scheme at \i = aT x using tadpole-improved one-loop matching 
factor [g], and then evolved to \i = 2 GeV using the four-loop RG-equation. 

Quark mass results are shown in Fig. |2| While 0(g 4 a) scaling violation should be present 
due to the use of one-loop matching factor, comparison of VWI and AWI masses for ud quarks 
suggests that these terms are small relative to the 0(a 2 ) term [Qj. Therefore we extrapolate quark 
masses linearly in a 2 . Possible effects of 0(g 4 a) terms are estimated from the ambiguity of the 
renormalization factor by either shifting the matching scale from jU = 1 /a to pi = %ja or using an 
alternative definition of coupling for tadpole improvement 

As already observed in Nf = 2 QCD values of the strange quark mass determined for 
either the K- or the <p -inputs, while different at finite lattice spacings, extrapolate to a common 
value in the continuum limit. Therefore the quark masses in the continuum limit is estimated from 
a combined fit to data with the K- and the <p -inputs. We finally obtain 

m^(n = 2 GeV) = 3.50(14)(+*j) MeV, m^(pL = 2 GeV) = 91.8(3.9)(+^) MeV. (4.2) 
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Figure 2: Continuum extrapolations of the up, down and strange quark masses obtained with the /T-input. 
For comparison, results for quenched and Nf = 2 QCD |Q| are overlaid. 

Dynamical up and down quarks reduce significantly the quark masses [0]. The effect of strange 
quark is less dramatic, and we do not see deviations from the Nf = 2 results beyond statistical 
errors. 

4.3 PS decay constants 

PS meson decay constants are estimated using matching factor determined by tadpole-improved 
one-loop perturbation theory. The results with K-input are 



f„ = 143.4(8.8) GeV, f K = 163.7(8.6) GeV, /*//„ = 1.140(17) 



(4.3) 



We recall that in our Nf = 2 QCD calculation, the magnitude of scaling violation was so large that 
we were not able to estimate values in the continuum limit [|p. The situation is much better in 
the present case and f K and fg turn out to be almost consistent with experiment. The errors are 
large, however. Furthermore, the ratio fa/ fx differs significantly from experiment. A long chiral 
extrapolation is a possible cause of the discrepancy. 

4.4 Sommer scale 

The Sommer scale in the continuum limit for A'-input reads 

r = 0.516(21) fm, (4.4) 

which is closer to a phenomenological value of 0.5 fm than the estimate ro ~ 0.541(17) fm in the 
N f = 2 QCD [§. 



5. Chiral fit using the Wilson #PT 

Our long chiral extrapolation is a dangerous source of systematic error. As a supplementary 
analysis, we try to fit meson masses using the #PT and compare it with the polynomial fit above. 

We use the #PT modified for the Wilson quark action (W^PT) JT0|], in which finite lattice 
spacing corrections to infrared chiral logarithms for the Wilson quark action are incorporated. The 
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Figure 3: Comparison of the W^PT fit and the polynomial chiral fits at j3 = 1.83. mf WI and the lattice 
spacing necessary to draw curves are set by polynomial chiral fits and the /T-input. 



Nf = 2 + 1 QCD W^PT formulae for the 0(a) improved theory has been obtained [11] up to the 
NLO. For our fit, we rewrite the formulae in terms of the AWI quark mass; 



ml = (x + Pgy) 



1 



1 + ^2 <J £ L ¥ (^i)A ¥ -H-l2L 46 (^)x-4L 5& (^)(x + P e y) 



1 



mo + te+Xyy-—— £ H*( g lg 2 2 )(x + P ¥ yf 2 , 



, (5-1) 
(5.2) 



where 6 = n,K and co = p,K*. In these equations Pg,A v and Hy(gf,g2) are geometric factors, / 
is the pion decay constant in the chiral limit of the LO formula, and 



x=^^(2m ud + m s ), y = ^(m u d-m s ), 
L 46 (At) = 2L 4 (aO +4L 6 (jU), L 58 (m) = 2L 5 (/x) +4L & (^i), 



(5.3) 
(5.4) 

16, 2 " ix 1 • (5 " 5) 
where Bo, L4, L5, L(, and Lg are the low energy constants in the continuum #PT. We note that up to 
NLO the vector formula has no lattice artifact whereas the PS formula has the term H which is of 

0(a 2 ). 

We find that while fit parameters vary largely among jackknife samples, and hence are not well 
determined, fitting curves themselves are stable. Also the W^PT fitting curves do not exhibit any 
significant difference from the polynomial fit, as shown in Fig. ^} Therefore, conclusions for the 
spectrum and quark masses derived from polynomial fits are not altered even if we use the W^PT 
fit. We suppose that a relatively large value of the up and down quark masses in our simulation 
makes it difficult to observe chiral logarithm behavior from data. 



6. Conclusions and future plans 

We have reported the results for the meson spectrum, light quark masses and other physical 
quantities obtained from our simulations in Nf = 2 + 1 QCD with Wilson-clover action. We find 
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that the meson spectrum is consistent with experiment, and the quark masses are smaller compared 
to often quoted phenomenological values. 

We regard the present result as a first step toward a fully satisfactory Nj ■ = 2 + 1 simulation. A 
major point to improve is the control of systematic errors due to quark masses still large compared 
to those in Nature. Larger volumes will also be needed for baryons, and even for mesons as quark 
masses are reduced. We hope to overcome these problems with the PACS-CS project [|l2|], using 



the improved algorithm provided by the domain decomposition idea [ ]13[ ] and the cluster computer 
PACS-CS developed at University of Tsukuba. 
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